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Abstract 



A direct insertion scheme for assimilating coastal acoustic tomo- 
graphic (CAT) vertical temperature sections into a multilevel, coastal 
primitive equation model for the Gulf of Sirte (Lybia) is investigated 
using computer simulation experiments. Although the model was 
developed for the whole Mediterranean Sea (MED), only a model 
sub-domain covering the Gulf of Sirte was used in this study. The 
model has realistic coastlines and bottom topography, and a coastal- 
following, curvilinear, nearly orthogonal, horizontal coordinate system 
with a horizontal resolution of about 10 by 10 km . The grid of the 
model was designed using a grid generation/focusing technique. The 
model has complete thermodynamics, second order turbulence closure, 
and 16 bottom-following (sigma) vertical levels. To generate the ’’true 
ocean” for this study, the model was first spun up for 30 days with 
the Levitus temperature and salinity and ECMWF wind climatolo- 
gies, and then run for one year and more. The last 60 days of this 
control run were taken to represent ’’actuality.” 

A series of assimilation experiments was carried out in which CAT 
temperature slices synthesized from different CAT configurations based 
on the ’’true ocean” were inserted into the model at various time steps 
to examine the convergence of this direct insertion scheme. In all the 
assimilation experiments, after a spinup for 30 days, the model was 
then integrated with CAT slices inserted daily to produce ’’nowcast 
fields” for periods of 30 and 60 days. 

Our results strongly indicate that the coastal model and the direct 
insertion scheme for CAT slices work well together. As time pro- 
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gresses, the CAT slices spread information out in nearly all directions 
in the Gulf of Sirte. The inserted data act as an initial convergence 
maker which forces the model to ’’lose” memory of the initial state 
and gradually converge to the ’’true ocean”. The global rms error of 
the nowcast temperature fields in the sub- domain considered in an 
experiment involving five CAT sections decreases by about 50% after 
30 days and 66% after 60 days of data assimilation. The rms errors 
of two other experiments involving only a single slice are reduced by 
33% after 30 days. In all three experiments, the rate of convergence, 
as measured by the rms error, is approximately constant. It is ex- 
pected that the convergence to the ’’true ocean” will be quicker with 
more frequent insertions and with more CAT sections. 

It is seen from these experiments that the temperature information 
spreads out faster over the shallow-water areas. This may be due 
to stronger mixing, horizontal shear and interaction with topography 
resulting in an increase of horizontal diffusion. The temperatures at 
the mixed-layer level converge more rapidly to the ’’true temperatures” 
in comparison to the surface-level temperatures. All the CAT slices 
have almost the same effectiveness in improving the nowcast fields 
regardless of their horizontal orientations. Our experiments also show 
that convergence is not linearly proportional to the number of CAT 
slices assimilated. Our simulation results confirm that the internal 
forcing provided by vertical slices derived by CAT is not localized in 
physical space. They induce immediate correlation between the fully 
nonlinear primitive equation model and the ’’true ocean”. 
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1 INTRODUCTION 



In the past ocean modeling tended to proceed separately from observations, 
mainly due to the lack of synoptic observations to properly constrain the 
model physics. Two new technologies in today’s Oceanography brings new 
large synoptic datasets. Satellites provide global maps of sea surface height 
(SSH), sea surface temperature (SST) and wind stress. The latter field is the 
most important surface forcing of the ocean circulation. The satellites give 
information on the physical characteristics of the ocean surface and the near- 
surface layer. On the other hand, Ocean acoustic tomography can provide 
quality 4-D data measuring the interior ocean structure. This underwater 
acoustic inverse technique uses travel time changes of sound pulses to map 
sound speed/temperature perturbation along various paths. The satellites 
and ocean tomography are highly complimentary ocean observing systems. 
The ability to efficiently assimilate these datasets into ocean models is key 
to ocean nowcasting and forecasting. 

Ocean tomography has traditionally been used in deep water for mapping 
the ’’ocean weather.” The adaptation of this technique from deep to shallow 
water poses many scientific and engineering challenges. New sensor systems, 
advanced sound propagation models, and optimal signal processing and map- 
ping techniques that will work in complex coastal environments must be de- 
veloped and tested at sea. Research into and development of coastal tomog- 
raphy is the main thrust of a collaboration between scientists and engineers 
at the Naval Postgraduate School (NPS) and Woods Hole Oceanographic 
Institution (WHOI). In August 1992, a first-generation coastal acoustic to- 
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mography (CAT) system was tested in the Barents Sea (Chiu et al., 1993, 
and Miller et al., 1993). This system consists of an electronically controlled 
sound source and a telemetered vertical hydrophone array (Von Der Heydt 
et al., 1992). Using such system, the two groups demonstrated that CAT is 
practical and can produce vertical slices of the coastal ocean temperature and 
sound speed fields at high temporal and spatial resolutions. An assessment of 
the potential performance of a coastal nowcasting system which assimilates 
CAT vertical slices into a state-of-the-art coastal ocean model is the focal 
point of a computer simulation study reported here. 

In the assimilation of oceanographic datasets, numerical ocean models 
play a critical role in the advection and propagation of the information con- 
tained in the data. Two outstanding questions were raised by Ghil and 
Malanotte-Rizzoli (1991): (1) Can information provided only at the sea sur- 
face be transferred dynamically into the deep oceanic layers, thus reconstruct- 
ing the deep circulation? (2) Can information provided only locally, in limited 
oceanic regions, be transferred to ocean areas far away from the data-dense 
region and which time and space scales are better estimated through the as- 
similation of local data? These issues, particularly pertaining to deep-ocean 
circulations, have been addressed by Ghil and Malanotte-Rizzoli (1991) as 
well as other experts working in the subject area. We will give a brief review 
of some of these previous computer simulation works which are somewhat 
related to our present work in studying the direct insertion approach and 
advection information process. The central difference between this work and 
the previous works is that this work focuses on nowcasting in the coastal 
oceans using CAT and a state-of-the-art coastal ocean model. 
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The problem of propagation of surface information to the ocean inte- 
rior was first studied by Hurlburt (1986) and Thompson (1986). They both 
inserted the altimeter data directly at every grid points into a two-layer prim- 
itive equation models for the Gulf of Mexico. Focusing on the dynamic trans- 
fer of information from the surface to the deep layer, Hurlburt (1986) demon- 
strated the success of the numerical ocean model used in the reconstruction 
of deep circulation. Hurlburt showed that approximately two updates per 
eddy cycle ( 57 days) are required for the assimilation to be successful and 
provide convergence to the reference ocean and the convergence is clearly 
accelerated further when updating every 20 days. Thompson (1986) focused 
on the geoid error as it affects the assimilation of altimeter data for mesoscale 
ocean prediction. He showed that errors in the geoid on spatial scales com- 
parable to the model gird resolution did not seriously degrade the forecast, 
even in dynamically active regions with large gradients in the geoid height 
where instability process has been observed to occur. Kindle (1986) used a 
one-layer reduced-gravity primitive equation model of the Gulf of Mexico to 
examine the methodology of incorporating satellite altimeter data into ocean 
forecasting model. His main result was that for a stationary circular eddy, 
approximately two track (either ascending or descending) across the eddy 
are sufficient to ensure adequate spatial resolution and an irregularly shaped 
eddy may require three or four tracks. 

In a two-part study, Malanotte-Rizzoli and Holland (1986 and 1988; here- 
after MRH) examined the transfer of information by advection from data lo- 
calized in space to other regions of the ocean using an idealized, mid-latitude, 
multilayer, quasi-geostrophic (QG), eddy-resolving model and a direct inser- 
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tion technique, weighting observations by their distance from the grid point 
being updated. In the first part (MRH, 1986) they used a quasi-linear, steady 
state model ocean and the assimilated data were density measured along vari- 
ous hydrographic sections. The main result of the study is that a local section 
can be quite effective in determining flow in regions far away from the data 
section if the flow is simple, steady, and quasi-linear. The most effective data 
sections are meridional, long and far away from the ocean’s western bound- 
ary. In the second part (MRH, 1988), they used a time-dependent, fully 
eddy- resolving model ocean and the assimilated data were sections with dif- 
ferent orientations. They found that, on the short time scale of mesoscale 
variability, all the data sections considered are equally ineffective. A single 
data section is quite effective in driving the model to the reference ocean if 
the data insertion process is carried out for time longer than the model equi- 
libration time. A single section of assimilated data is completely ineffective 
in driving the model toward reality over time scales comparable with the 
mesoscale variability. Their result also show that the model climatological 
mean after ten years of data assimilation becomes extremely similar to the 
climatological mean of the reference ocean. 

For this present study, we use a highly-nonlinear, primitive equation, 
eddy-resolving, coastal model for the Gulf of Sirte (the Gulf of Lybia). 
The model uses curvilinear coordinates with coastal-following and nearly- 
orthogonal grid. The model physics include fully active thermodynamics 
and a second order turbulence closure. This coastal model was developed by 
adapting the Princeton Ocean Model (POM) for the Gulf Stream (Mellor and 
Ezer, 1991) to the Mediterranean Sea (MED). Previously, Mellor and Ezer 
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(1991) used POM and an optimal-interpolation altimeter-data assimilation 
scheme to study the quality of the nowcast and forecast fields in the Gulf 
Stream region. Here, we use the model to study the potential of using CAT 
for coastal ocean nowcasting. Our assessment is focused on the simple direct 
insertion approach. An interesting question is that what configurations of 
CAT in a shallow water region such as the Gulf of Sirte (Gulf of Lybia), 
where both topographic and nonlinear effects are strong, are most effective, 
i.e., lead to a large and widespread improvement of the nowcast fields. Our 
emphasis here is on the understanding of horizontal propagation of the infor- 
mation obtained by CAT, the effect of data insertion intervals on the nowcast 
fields, and the convergence rate associated with the direct insertion scheme 
using this state-of-the-art curvilinear, multilevel primitive equation coastal 
ocean model for the Gulf of Sirte. 

The report is organized as follows. In the next section, we give a descrip- 
tion of coastal acoustic tomography (CAT) with a focus on the inversion of 
tomographic data to obtain vertical temperature slices. The coastal ocean 
circulation model for the Gulf of Sirte is described in Section 3. The model 
was first used to simulate the ’’true ocean” and, subsequently, it was used to 
produce nowcast temperature fields with the assimilation of synthetic CAT 
slices for various CAT configurations. The direct insertion technique is also 
discussed in Section 3. In Section 4, the numerical experiments and the 
results are discussed. Section 5 gives a summary of our findings. 
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2 COASTAL ACOUSTIC TOMOGRAPHY 



In contrast to deep water propagation, the arrivals from the multipaths due 
to a pulse excitation in a shallow-water environment tend to overlap in time. 
It is therefore difficult to resolve them individually with omni-directional 
or almost omni-directional receivers which are typically used in deep water 
tomographic applications. Vertical line arrays and beamforming techniques 
can be used to overcome this difficulty. The observed perturbations of the 
beamformed individual ray/mode travel times constitute the database of the 
coastal tomographic inverse scheme developed by Chiu et al. (1993). Their 
scheme is discussed next in detail. 



2.1 A Hybrid Ray/Mode Inverse Method 

Establishing the mathematical relations between the observations and the 
unknown structure is the first step in the development of any inverse scheme. 
In our case, the observations are the changes in ray and mode travel times, 
Stj and 8t™, and the unknown ocean variable is perturbation of sound speed, 
8c, i.e., the deviation from a reference sound speed field c. The reference 
ocean is allowed to be range dependent in our formulation. 



2.1.1 Ray Travel Time 



The relation of the travel time of a pulse signal along a raypath, t r j , to the 
speed of sound c is well known. The relation is 





ds 



( 1 ) 
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where x = (z, j/, z) is the position vector and s is the arclength along the path 
trajectory Tj connecting a source and a receiver. Replacing c by c + 6c and t r - 
by tj + Stj and assuming the change in travel time due to path perturbation 
is of 2nd order, we obtain, through linearization, 




(2) 



2.1.2 Mode Travel Time 



The derivation of the relation between St ™ and Sc is less straight forward. The 
uniformly valid, asymptotic coupled-mode solution of Desaubies et al. (1986) 
shall be our starting point. This solution was obtained by expanding both 
the magnitude and phase of the modal horizontal structure in a perturbation 
series in powers of a small parameter p, where p is the ratio of the wavelength 
to the characteristic horizontal scale of the oceanic variability. The horizontal 
phase of the n-the mode, to 0(p 2 ), is given as 



a f T < \N 1 A m n ( r ?k-0 u 

0»=/ K»(r,w) 1 Tt \ 2 ? \i dr 

Jo 2 ^ K n 2 (r,u) - K m 2 (r,oj) 



( 3 ) 



where r is the range along the horizontal mode path connecting a source and 
a receiver, k„ and K m are the horizontal wavenumbers of modes n and m, 
respectively, and A mn is a coupling coefficient that can be equated as 



„ f°° 1 _ . A dZ m (z;r,u) v \ 

Jo p(r, z) or 



dZ n (z;r,u>) 



]dz (4) 



o ''i } J C\ 

or or 

In the depth integral of (4), p is density, and Z n and Z m are the nth and 
mth modes, respectively. 

Analogous to the ’’static ray approximation,” we assume that the varia- 
tions in the mode functions produce higher-order changes in the horizontal 
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phases comparing to those caused directly by the variations in the horizontal 
wavenumbers due to a perturbation in sound speed. This ” static mode ap- 
proximation” was introduced by Rajan et al. (1987) in their inverse method 
for obtaining geo parameters. Under this assumption, the variation of the 
modal phase integral (3) leads to 



86 , 



TOO 

n = / ( 8li n + fmn&Kn H" 9mn^ ^m)df 

^ 0 m rr> 



(5) 



with 



_ 1 (K n 2 + K m 2 )A mr ? 

Jmn — 2 { - Kn 2_ Rm 2y 



9 ran — 



A m n 



(6) 

(7) 



(«n 2 - «7n 2 ) 2 

where the modal quantities with an ’’overbar” are those associated with the 
reference ocean. From a stationary phase consideration, modal travel time 
and phase are related by 



dd n 



t m = 
n du 



( 8 ) 



It follows from (8) and (5) that the variation in modal travel time is 

386 ' 



= a. 



[ r 8&t C n U Jmn f ^ 

= /o ( ^7 + ?-a7 ii ' c “ + ? / - 



df„ 



88 K n 

'~KT 



(9) 



8(Jmn i~ 88n r 



+E^ + Ew„ 

™ 8u „ du 



)dr 



To find the relation between 8t ™ and 8c , we need the relations of 8tc n 

r\ r 

and hence to 8c. Using the governing vertical-structure equation and 
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orthonormal condition for the normal modes, Raj an et ah (1987) showed 
that the perturbation of AC n due to a small Sc is 

TOO 

6k u = / q n 6cdz 
Jo 




with 

— Cc> 2 1-2 

Note that the ’’static mode approximation” was also imposed in obtaining 
the eigenvalue perturbation equation. Finally, a substitution of (10) in (9) 



gives 



K 



- /:/: < 



°° r dq n (r,z,u>) 



du> “ du 

+g m n(r, z,uj)q m (r, z,uj)]}6c(x)dzdr 



+ H Q-[/mn(r, Z, w)9„(r, z, u) 



( 12 ) 



Previously, Shang (1989) has developed a tomographic inverse technique 
based on adiabatic mode theory. His formula for modal travel time per- 
turbation has only the first term in the integrand. Here, we have included 
the effects of mode coupling in the travel time change. 



2.1.3 Optimal Inverse Solution 



We represent Sc in the vertical slice traversed by the rays and modes as a 
Fourier series, i.e., 



with 



1/2 K/2 

Sc = E E a.ijk exp[\^-T(zX 0 x + kV 0 z)\ 

i=-I/2 k=—K/2 



r 

K °- X 



V 0 = — 



27T 

~Z 



(13) 

(14) 

(15) 
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where X and Z are the horizontal and vertical extents of the slice, respec- 
tively, and (/ + 1) x (K + 1) is the total number of spectral components 
used. This truncated Fourier decomposition of Sc is equivalent to gridding 
(i.e., discretizing) the physical space (a:, z ) into (/+ 1) x (K + 1) points with 
a horizontal grid spacing of ^ and a vertical grid spacing of Using the 
Fourier decomposition, the inverse problem can be cast into simple linear 
algebra. 

Putting the Fourier coefficients a, Vs in a column vector, a, and collecting 
all the measured ray travel time changes along resolved multipaths connect- 
ing sources and receivers in a column vector, 8t r , the corresponding data- 
unknown relation becomes 



St r = G r a + e r (16) 

where G r is the ray data kernel matrix, calculated by numerically integrating 
(2) for each of the Fourier components. Similarly, the data-unknown relation 
corresponding to the modes can be cast as 

8t m = G m a + e m (17) 

where the data vector Si™ contains the measured mode travel time changes 
and the mode data kernel matrix G ra is computed from (12). In general, the 
data are contaminated by experimental noise, a combination of measurement 
and model errors. We represent the contamination in (16) and (17) by the 
noise vectors e r and e m , respectively. 

There are several ways to derive the linear, minimum mean-square-error 
estimate, i.e., inverse solution (Cornuelle, 1985, and Chiu, 1985). Assuming 
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that a, e r and e m have zero means and are uncorrelated, it can be shown that 
the minimum mean-square-error solution can be obtained by minimizing the 
following objective (or cost) function: 

5(a) = - G r a) T C£ -1 (£f - G r a) 

+!(&”* ~ G m a) T C?-\8t m - G m a) 

+— a 7 'C a l a (18) 

where C£, C™, and C a are the covariance matrices of e r , e m and a, respec- 
tively. Furthermore, the error covariance of the solution is simply the inverse 
of the Hessian (the matrix of second derivatives) of S (Chiu, 1985). One 
can see that the objective function is constructed by a combination of two 
data constraints (associated with rays and modes, respectively) and a statis- 
tical information constraint. It is worth mentioning that if e r , e m and a had 
normal statistics, the minimum of (18) would coincide with the mode of the 
a posteriori distribution, which is the probability distribution function of a, 
given 6t r and 6t m . In such case, the minimum mean-square-error estimate 
and the mode of the a posteriori distribution estimate would be identical 
(Chiu, 1985). 

The minimum mean-square-error inverse solution, a, and its error covari- 
ance, C c , follow from the evaluations of the first and the second derivatives 
of S with respective to a. They are 

a = C £ (G rT C^- 1 M r + G mT Cf" 1 6i m ) (19) 

C £ = (G rT Cr 1 G r + G mT Cf- 1 G m + C a “ 1 ) -1 (20) 

It is straight forward to go back to the physical space ( x,z ) from the 
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Fourier domain. The inverse solution for sound speed perturbation, 8c, and 
its error covariance are obtainable through applications of multi-dimensional 
discrete Inverse Fast Fourier Transforms (IFFT) to their Fourier-space coun- 
terparts. Ocean temperature perturbation is linearly proportional to sound 
speed perturbation. Thus, the calculation of temperature estimates from 
sound speed estimates is straight forward. 

3 COASTAL OCEAN CIRCULATION MODEL FOR THE GULF 
OF SIRTE 

3.1 The Model 

A coastal primitive equation model with a fully active thermodynamics is de- 
veloped for whole Mediterranean Sea (MED) based on the Princeton Ocean 
Model (POM) for the Gulf Stream (Mellor and Ezer, 1991). This is a 
three-dimensional MED model with a horizontal coastal-following, curvilin- 
ear nearly orthogonal and bottom-following, sigma, vertical coordinate sys- 
tem. The horizontal coastal-following, curvilinear nearly orthogonal model 
grid is designed using the EAGLEView software package developed at Missis- 
sippi State University (MSU). While the package allows for a high-resolution 
handling of the open boundary, coastlines and Gulf regions, a reduction of 
computational expense is still achieved (Ly et al., 1993). The model has 
free surface dynamics and it uses a second order turbulence closure scheme 
for the mixed layer and a split mode time step. The model has one open 
boundary in the Atlantic ocean about 200 km west of the Strait of Gibral- 
tar. Only the model sub-domain covering the Gulf of Sirte is used in this 
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data assimilation study. The average horizontal resolutions for the studied 
region is approximately 10 km. 

The model equations governing the velocity, surface elevation, salinity, 
and temperature fields in the ocean have been described by Blumberg and 
Mellor (1983; 1987). The equations are written in a system of Cartesian 
coordinates with x eastward, y northward, and 2 upward. The free surface is 
located at z = rj(x,y,t) and the bottom is at z = — H(x,y). The equations 
are 



!E + u.VU + W^- + 2fix[/ = — —VP + + F (21) 

at dz po dz dz 



dP_ 

dz 



= -99 



( 22 ) 






(23) 



ddi - dd; d ddi 

+u . vtl+w . )+F < 



(24) 



The density is computed using the equation of state (Mellor, 1991) in the 
general form 

P = p${S,0,p) (25) 

The Coriolis force is denoted as 2Q x U, where Q is the earth rotation 
vector, U is the horizontal velocity vector with components ( U , V), V is the 
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horizontal gradient operator, p Q is the reference density, p is the in situ den- 
sity, g is the gravitational acceleration, P is the pressure, K m and Kh are 
the vertical turbulent exchange coefficients for momentum, and for heat and 
salt, respectively. In equation (24), 9{ may represent mean potential temper- 
ature, 0 (or in situ temperature for shallow water application) or salinity, S. 
The potential density used here is an approximation, since it excludes the 
effects of pressure variations (Blumberg and Mellor, 1981). Here, F(F x ,F y ) 
in equation (21) and F$ t in equation (24) are the horizontal mixing terms 
(see Blumberg and Mellor, 1987). 

The momentum conservation and diffusion equations (21) and (24) con- 
tain the vertical turbulent exchange coefficients which are determined by a 
second order turbulence closure scheme of Mellor and Yamada (1982). The 
turbulence scheme is characterized by equations for turbulent kinetic energy 
(TKE), q 2 / 2, and for the turbulent mixing length, i. The two equations can 
be written in the same form (see Ly, 1992a), so that Qi is either q 2 / 2 for 
TKE or q 2 £ for the turbulent mixing length in the equation: 



dQi +U-VQi+W dQi 



d . dQi 



,dU 



dV , 



9 dp 



dt 






dz dz 



dz 



dz 



po dz 



(26) 



B i 



where 



K mQ = 2K m , I< hQ = 2 AT, K q = 2/£, F Q = F q (27) 

when Qi represents (q 2 / 2) 
and 



KmQ = £E\K m , KhQ = £E\Kh-> K Q = W , Fq = Ft 



17 



(28) 



when Qi represents (q 2 £). Also in the above equations, W is a ”wall proxim- 
ity” function ( Mellor and Yamada, 1982) and B\ is an empirical constant. 

The first term on the right-hand side of (26) is the diffusion contribution, 
the next two terms represent shear production, the next is the buoyancy term, 
the next term represents dissipation, and the last represents the horizontal 
mixing. The vertical turbulent exchange coefficients (TEC) K m , Kk , and K q 
are defined as K m = £qS m , Kh — ZqSh, K q = ZqS q where Sm^Sh, and S q 
are the stability functions (see Mellor, 1973; Mellor and Yamada, 1982). 

At the free surface, z — rj(x, y), the surface wind stress, heat, and salinity 
fluxes are prescribed. At the bottom z = — 77(x,y), zero heat and salinity 
fluxes are used. Turbulent mixing length, £, vanishes at the bottom. At 
land boundaries the condition of no diffusive fluxes of any property across 
the interface is used. Sigma coordinates can adequately model domains like 
the Gulf of Sirte with large bathymetric irregularities. Mode splitting of the 
barotropic and baroclinic modes are imposed in the model to save computer 
resources. 

3,2 Model Parameters and Initialization 

The Gulf of Sirte off Lybia is located in the southern middle part of the 
Mediterranean (MED). The region of interest in this study is bounded by 
latitudes 30 and 35 N and longitudes 15 and 22 E. The eastern part of the 
Gulf of Sirte is a deep region where the 100 — m isobath is less than 20 
km from the coastline. The southern and western parts of the Gulf are 
considerably shallower. The bottom depths along the CAT sections used for 
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the assimilation are typically from 200 m to 2000 m. 

In the model, the Gulf of Sirte and its surroundings have a horizontal res- 
olution of about 10 x 10 km. There are 16 bottom-following (sigma) vertical 
levels with increased resolution in the mixed layer and lower resolution in the 
deeper layers. The external (barotropic) mode time step is 15 s, and the in- 
ternal (baroclinic) one is 10 minutes, satisfying the Courant-Friedrichs-Levy 
(CFL) computational stability criterion. The model bathymetry is obtained 
from the global DBDB5 bathymetry dataset with 5' x 5' resolution (National 
Geophysical Data Center, 1985) and is interpolated to the model grid. The 
bathymetry is filtered with a Shapiro (1970) filter to remove high frequency 
noise. The model is initialized with the Levitus temperature and salinity 
climatology (Levitus, 1982), and wind stress from the ECMWF (Ly et al . , 
1992b) climatology. In all our computer simulation experiments, the model is 
spun up for 30 days diagnostically (temperatures and salinities held constant) 
before making the prognostic run. 

3.3 Data Direct Insertion Technique 

There axe various data assimilation approaches. The simplest and the most 
computationally tractable is the direct insertion approach which is used in 
our investigation. In the following discussion of the direct insertion technique, 
we use superscripts /, o, and a to denote the forecast, observed, and analyzed 
(the best first guess) ocean fields, respectively. We also assume that CAT 
data is available at time intervals of At. 

Generally, the ocean field predicted at time t , prior to the assimilation of 
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the data obtained at t, is related to the analyzed ocean field at time t — At 

by 

T[ = r(Tf_ At ) ( 29 ) 

where T is the mathematical operator describing the ocean dynamics, which 
in our cast is nonlinear. In the direct insertion method, the inserted data can 
be viewed as a ’’perturbation” introduced at a specific location {x s ,y s , z 3 }. 
We insert data at all time stepes. The vector T represents the discretized 
ocean variables such as velocity and temperature at the ocean model grid 
points. An observational vector can be thought of as consisting of a predicted 
and a perturbation part e as: 



T? = T[ + e (30) 

Then in the assimilation experiments we have 

r(T?_ it ) = T? = T{ + £ (31) 

at specific locations {a;,y,z} = {#«,?/«, 2 S }. And 

r(Tf.At) = Tf (32) 

at all other locations. Where £ is a difference between the predicted and 
observed value (a perterbation) and grows in time and space. It propagates 
and diffuses away from the CAT slice section. It is noted that the ’’true 
ocean” can be produced by replacing the predicted value T[ with an observed 
T° at spatial and temporal resolutions which are dense and frequent enough 
to force the difference e toward zero. 



20 



The analyzed field at time t — At in the insertion experiments can be 
constructed by combining the forecast field and the observed field at t — At, 
that is 

TLa, = (I - W)T{_ At + wt;. 41 (33) 

where I is an identity matrix and W is a diagonal weighting matrix. The 
diagonal elements of W consist of zeros and ones. The element is one if the 
corresponding model grid point is traversed by the CAT sections and if the 
corresponding variable is temperature. The element is zero otherwise. The 
construction of the analyzed field using equation (33) and the prediction of 
the forecast field using equation (29) is repeated at every data assimilation 
time step. 

4 NUMERICAL EXPERIMENTS AND DISCUSSIONS 

In addition to the model run that generated the ’’true ocean” (which will 
be referred as the TR experiment), four other model runs were performed. 
The latter runs will be referred as the EX n experiments with n = 0,1, 2, 3. 
EXo is an experiment with no data inserted. This experiment provides the 
baseline for judging improvements when data are inserted. In EX i (fig- 
ure 10), the model was run with a single CAT section inserted daily. This 
CAT section has a northeast-southwest (NE-SW) orientation. In EX 2 (figure 
17), a single CAT section was used again but, this time, the section has a 
northwest-southeast (NW-SE) orientation. In EX 3 , the model was run with 
the insertion of a total of five CAT sections. The five sections include the 
previous two diagonal slices as well as three other slices enclosing a coastal 
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volume making a CAT parallelogram (figure 24). These experiments are 
discussed in detail next. 

4.1 The ’’True Ocean” Experiment TR 

The scheme of the twin experiments in this study is shown in Table 1. The 
central issue is the convergence problem in using a fully nonlinear primitive 
equation coastal model with the direct data insertion technique. The global 
volume rms errors between the temperature fields of various experiments and 
the ’’true ocean” can be measures of the convergence in the domain studied. 

The model is spun up for 30 days and then run for a year with ECMWF 
winds to produce the ’’true ocean” TR. One year is run to make sure the 
model ’’loses” the initial state (the loss of memory of the initial state). Fig- 
ures 1, 2, and 3 show the ’’true ocean” temperature fields for day one at the 
surface (level 1), the mixed layer level (level 6), and the bottom level (level 
15), respectively. Figures 4-6 show the ’’true ocean” temperature fields at the 
above levels for day 30. Hereafter, in the TR run, days 1 and 30 correspond 
to days 361 and 390, respectively. After one year of simulation the model is 
stable with the annual ECMWF winds. From figures 1 and 4 it is seen the 
SST fields are very slightly developing after 30 days (after one year of model 
simulation). The SST field in the shallow waters (the southern, southwest- 
ern, and southeastern regions of the Gulf) stay very much the same after 30 
days of simulation. Temperatures in these regions are about 13 deg C . In the 
central and northern parts of the region, the SST field is slightly different. 
The region of 12 deg C has developed towards the northwest after 30 days 
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simulation. The colder water of the deeper layer has moved up to the surface 
by the mixing process in the central and northern parts (the deepest region 
of the MED with depths of about 4000 m). 

Figures 2 and 5 show the temperature field of the ’’true ocean” at level 6 
for one and 30 days, respectively. It is noted that level 6 of sigma coordinates 
corresponds to depths from the subsurface to 100 m. The temperature fields 
in all parts of the region stay very much the same after 30 days simulation. 
They have changed little compared to the surface level (figs 1 and 4). The 
temperature fields are almost identical at the bottom level (level 15) after 30 
days simulation which are shown in Figures 3 and 6. 

As expected, the 30-day run changes slightly the surface temperature 
field. In the deeper layers (levels 6 and 15), the temperature fields stay 
almost the same after the 30-day run. 

4.2 The No Data Assimilated Experiment EX 0 

The model is spun up for 30 days with Levitus temperature and salinity and 
ECMWF wind climatologies. The no data assimilated experiment EXo is 
performed by running the model for 30 days without the CAT data insertion. 
The temperature fields without data insertion at days one and 30 are shown 
in Figures 7 and 8, respectively. The temperature is developed after the 30- 
day run (Figure 8), but still close to the day one temperature when no data 
is inserted. 

The effectiveness of the insertion of CAT data can be illustrated in Fig- 
ure 9 which shows the global volume rms values of the differences between 
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this experiment without data assimilation and the comprehensive experiment 
with data insertion EX 3 . Figure 9 shows that at an initial time (day 1) the 
differences are very small. The differences are close to zero at day 1 because 
the temperature fields of EXo and EX S are the same at day zero. The fields 
are more different after more days of insertion. The growth of the global 
rms value of the differences must be understood as a ’’good sign” of the data 
insertion process. Here, the inserted data acts as an initial convergence fac- 
tor to force the model to "lose” the initial state quickly and converge to the 
’’true ocean” (see the experiments EX j , EX Z and EXs). The global rms dif- 
ferences grow very fast in the first ten days, slow down in the next ten days, 
and becomes stable at about day 30. A measure of the time scale for loss of 
predictability of the model is T ime pre di c tabiiit y =~ 30 days in the experiment 
EX 3 and other experiments (not shown). 

4.3 The NE-SW Assimilation Experiment EX t 

The NE-SW CAT slice section can be seen in Figure 10 by a chain of 11 
squares which has a distance of 300-350 km. The model is spun up for 30 
days with temperature, salinity, and wind climatologies which are described 
above. Then the EX t experiment is performed by the insertion of CAT 
data for 30 days in the NE-SW slice at one day time step and at all levels. 
Figures 10-12 and 13-15 show the EX t temperature fields at the surface (level 
1), mixed layer level (level 6), and the bottom level (level 15) for days one 
and 30, respectively. Comparing Figure 1 of the ’’true ocean” and Figure 10 
of the NE-SW slice assimilated at the surface after one day of insertion, it 
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is seen that the two surface temperature fields are very different from each 
other, both in shallow and deep water. Here, the model accepts the CAT 
NE-SW data slice and spreads it out in both directions but more towards the 
east. Figure 13 shows the field after 30 days of assimilation. Figure 13 also 
shows the spreading is increased in the shallow water which may result from 
stronger mixing and horizontal shear and its interaction with topography. 
The stronger horizontal shear increases the horizontal diffusion [terms F and 
F$ ; in equations (21) and (24)]. This is made clearer by comparing Figure 
14 of day 30 assimilation at level 6 with Figure 2 of the ’’true ocean” at 
the same level. Here, the CAT data from the NE-SW assimilated slice is 
accepted by the model and is spread out towards the southeast where the 
Gulf has shallower average depths. In the bottom layer the Figure 6 of the 
’’true ocean” and the Figure 15 of the NE-SW assimilation are very much 
the same. 

The convergence of the NE-SW CAT slice assimilation experiment to- 
wards the ’’true ocean” can be measured by the global rms-error over a 
volume of the entire domain studied. Figure 16 shows the rms-error between 
the ’’true ocean” and the NE-SW slice experiment EX} for 30 days of as- 
similation. The rms-errors are 0.95, 0.925, and 0.64 for days 0, 1, and 30, 
respectively. This means that the rms-errors of the NE-SW slice assimilation 
decrease 3 and 33 % after one day and 30 days of assimilation. The model 
and assimilation scheme work very well in this case. The rms-error curve is 
almost linear and it is expected that the convergence to the ’’true ocean” 
would be even better if the insertion were made more often and at higher 
horizontal and vertical spatial resolution. 



25 



4.4 The NW-SE Assimilation Experiment EX 2 



This experiment is very similar to the NE-SW CAT slice above. The NW-SE 
CAT slice can be seen in Figure 17 by a chain of 11 squares over a distance 
of 300-350 km. In this experiment, the model is spun up for 30 days with 
the Levitus temperature and salinity, and the ECMWF wind climatologies. 
Then the model is run for 30 days with the CAT data insertion along the 
NW-SE slice at one day time interval for all levels (16 level) to perform the 
assimilation experiment EX 2 . 

Figures 17-19 and 20-22 show the EX 2 temperature fields at the level 1 
(surface), level 6 (mixed layer), and level 15 (bottom layer) for days 1 and 
30, respectively. In general, it is seen from these figures that temperature 
fields are developing towards the ’’true ocean”. As in the experiment EXj , 
comparing Figure 1 of the ”true ocean” and Figure 17 of the NW-SE slice 
assimilated at the surface, shows that the two surface temperature fields 
are very different from each other everywhere in the region. The surface 
temperature field after 30 days of CAT assimilation is shown in Figure 20. It 
is seen that the model accepts the CAT slice data and spreads it out along 
both sides of the NW-SE section at the surface. It is clear from Figure 5 
of the ’’true ocean”, Figure 18 of day one, and Figure 21 of day 30 that 
at level 6 the temperature of the NW-SE slice is shifted and more spread 
out northwards in the mixed layer. As mentioned in the EX t , this may be 
the result of the strong horizontal shear in the NW-SE slice region which 
increases the horizontal diffusion [terms F and Fe x in equations (21) and 
(24)]. In the bottom layer the temperature fields of the ’’true ocean” at day 
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30 in Figure 6 and of the NW-SE slice assimilation experiment EXg at day 
30 in Figure 22 are very much the same. 

The convergence of the temperature field to the ’’true ocean” in the NW- 
SW CAT slice experiment EXg can be measured by rms-errors in Figure 23. 
This figure shows the global rms-error over the volume of the entire domain 
between the ’’true ocean” and the the NW-SE slice experiment EXg. It is in- 
teresting to note that the global rms-errors of EX t and EXg have nearly the 
same distribution for 30 days of assimilation which are 0.95, 0.925, and 0.065 
for days 0, 1, and 30, respectively. As in EXj, the rms-errors of the NW-SE 
experiment decrease 3 and about 33 % after one day and 30 days of assimila- 
tion, respectively. One of the reasons for the similarity in the two rms-error 
distributions is that two slices (NE-SW and NW-SE) have the same length 
(300-350 km) and are very symmetric. As in the experiment EXj , the model 
and the assimilation scheme work very well together in this case. The model 
is greatly improved (33 %) after 30 days of CAT insertion. It is also seen that 
the behavior of the global rms-errors is nearly insensitive to the change in the 
CAT data section location in energetic regions of the Gulf in the experiments 
EX i and EXg. These two slices have almost the same effectiveness in improv- 
ing the model errors on the global space scale. These experiments also show 
that the internal forcing provided by the single slice of inserted CAT data is 
not localized in physical space in inducing immediate correlation between the 
fully nonlinear primitive equation model and the ’’true ocean”. This result 
is very different from the results of Malanotte-Rizzoli and Holland (MRH, 
1988) who found that the single section of inserted data is too localized in 
physical space to induce any immediate correlation between a QG model and 
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the reference ocean. The reason for these different results is the differences 
not only in data and models used, but also in our smaller coastal domain. 
The rms-error curve is almost linear in this case and it is expected that the 
convergence to the ’’true ocean” will be better if the temporal and spatial 
resolutions of the CAT insertion is higher. 

4.5 The Combined Slices Assimilation Experiment EX 3 

The more comprehensive experiment is EX 3 which is more realistic for coastal 
acoustic tomography. The combined slices in experiment EX 3 include the 
EXj slice and EX 2 slice together and three other CAT slices. All these slices 
together make a parallelogram without the southern side. The EX t and EX 2 
slices are two diagonals of this CAT parallelogram which is seen in Figure 24. 
The lengths of the western and eastern sides of the parallelogram (which will 
be called LON LEFT and LON RIGHT) are about 200 km, and the length 
of the northern side (which will be called LAT slice) is about 250 km . As in 
two previous experiments, the model is spun up for 30 days with the Levitus 
temperature, salinity, and ECMWF wind climatologies. Then the CAT data 
is inserted into the model, as described in the scheme (31) and (32), for 30 
days at one day time intervals and at 16 levels to perform the combined slices 
assimilation experiment EX 3 . 

The temperature fields of the combined CAT slices after one day of as- 
similation at the surface (level 1), the mixed layer level (level 6), and the 
bottom levels (level 15) are shown in Figures 24-26. Figures 27-29 show the 
temperature fields of the above levels for day 30. It is seen from these figures 
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that the temperature fields of all levels at day one are very ’’far” from the 
’’true ocean”. As in two other experiment, the temperature fields in this 
experiment of levels 1 and 5 are closer to each other than to the fields of the 
bottom layer. This is a result of the surface layer and the mixed layer levels 
being much closer with stronger mixing, especially in shallow water. After 
30 days of the combined slices assimilation, the temperature at the surface 
(Figure 27) is spread out in all directions, but stronger in the southeastern 
corner and in the southern part of the region where the water is shallower. 
The mixing, horizontal diffusion, nonlinearity, and interaction with topogra- 
phy play an important role in the advection process of temperatures. At day 
30 the SST field of this experiment (Figure 27) converges well to the ’’true 
ocean” of Figure 4. It is seen in all three experiments that the temperatures 
of the mixed layer level (Figures 14, 21 and 28) at day 30 converge very well 
to the ’’true ocean” in comparison with the surface levels. This results from 
the fact that the mixing is strongest inside the mixed layer where the vertical 
turbulence has the largest activity. The turbulence exchange coefficients K m 
and Kh in (21) and (24) have maximum values inside mixed layers and min- 
imum values at the surface (see Ly, 1991). Inside the mixed layer both ’’top 
down” and ’’bottom up” processes must be active while only the ’’bottom 
up” process is active at the surface. These processes are important in our 
cases because of the CAT data insertion by vertical slices. The temperature 
field at level 15 after 30 days of assimilation in Figure 29 is almost identical 
with the ’’true ocean” in Figure 6. 

The convergence temperatures can be seen in Figure 30 which shows the 
global rms-error over a volume of the entire domain between the ”true ocean” 
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and the the temperature in this experiment EX $ . As expected, the conver- 
gence of the temperature field to the ’’true ocean” is best in this experiment 
in comparison with the two experiments above. The rms-errors drop dra- 
matically in the first ten days of assimilation from 0.92 to 0.72 deg C. This 
makes the convergence to the ’’true ocean” of 22 % after the first ten days of 
assimilation. The rms-errors drop slower in the next 10 days from 0.72 at day 
10 to 0.59 at day 20. The convergence in this second ten days of assimilation 
is 18 %. In the last ten days of assimilation the rms-errors drop slowest from 
0.59 to 0.50 at day 30 which makes the convergence of 15 %. Totally, the 
global rms-errors drop is about 50 % after 30 days of data insertion. The 
model and the assimilation scheme work very well together. This shows that 
inserting the ’’right” data into a ’’wrong” model tends to drive the ’’wrong” 
model closer to the ’’truth” and to reduce global rms-errors. The inserted 
CAT data acts as an initial convergence maker in the model. This conver- 
gence grows very fast. After 30 days of insertion the convergence is increased 
by almost 50 %. This also shows that the ’’true ocean” and the model with 
CAT data insertion are strongly correlated. As in the two above experiments, 
the rms-error curve is very close to linear and it is expected the convergence 
to the ’’true ocean” will be better if the temporal and spatial resolution of 
the CAT insertion is higher. It is interesting to note that convergence are 
not proportional to the number of CAT slice sections in an experiment. The 
experiment EXi has a convergence of 33 %; EX 2 has a convergence about 
33 %, but EX S has about 50 %). 

It is interesting to see the convergence of this experiment EX$ to the 
’’true ocean” after another 30 days of data insertion. Figure 31 shows the 
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temperature field at the mixed layer level (level 6) at day 60 in comparison 
with the ”true ocean” at this day. The temperature converges significantly 
to the ’’true ocean” at day 60. The temperature fields at this level are almost 
identical. The rms-error distribution for days between 31 and 60 is shown 
in Figure 32. The rms-errors between the ’’true ocean” and the temperature 
field of this experiment are reduced by 17 % after 30 days (31-60) of insertion. 
Totally, after 60 days of CAT data insertion the rms-errors are reduced by 
66 %. The model is improved by 66 % after 60 days of assimilation. 

Three other assimilation experiments with CAT insertions along the north- 
ern (LAT), western (LON LEFT) and eastern (LON RIGHT) sides of the 
above parallelogram (Figure 24) have been carried out to study the effective- 
ness of various local slices. The model spreads these CAT slices out differ- 
ently, and this strongly depends on the horizontal diffusion, mixing process, 
and nonlinear interaction with topography. At the surface, after 30 days of 
assimilation the LAT slice is spread out along both sides, but more towards 
the north (not shown). The LON LEFT slice is spread out greater in the 
northern part towards the east, after 30 days (not shown). And at the sur- 
face the LON RIGHT slice, after 30 days, is spread out greater in the middle 
towards the west (not shown). Although these CAT slices have various tem- 
perature advection, but the their global rms-errors distributions are almost 
the same (not shown). Their rms-errors are reduced by 27 % after 30 days 
of assimilation. These experiments confirm that there is no most effective 
slice between different CAT slice sections. All slices have almost the same 
’’weight” in forcing the coast model in converging to the ’’true ocean”. 
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5 SUMMARY 



Coastal acoustic tomographic (CAT) data is used in ’’twin” experiments by 
direct data insertion. A series of assimilation experiments has been carried 
out and discussed in which data are inserted along various CAT slices into a 
coastal ocean primitive equation model for the Gulf of Sirte. The model is 
developed for the whole Mediterrainean (MED) Sea and is used for studying 
a direct insertion technique. The model has realistic coastlines and topogra- 
phy. This model also has a coastal-following, curvilinear, nearly orthogonal 
horizontal coordinate system, fully active thermodynamics, a second order 
turbulence closure, 16 bottom-following, sigma vertical levels, and an aver- 
age resolutions of about 10 by 10 km . The model was spun up for 30 days 
with Levitus temperature and salinity and ECMWF wind climatologies and 
then run for a year to generate the ”true ocean”. In the assimilation experi- 
ments, after a 30-day spinup, the model was run for 30 days with an insertion 
time interval of one day at all 16 vertical levels (except the comprehensive 
experiment EX 3 which was run for 60 days). 

The results indicate that the model works well with the assimilation 
scheme to accept data slices and spread it out nearly in all directions in 
the Gulf of Sirte. The inserted data acts as an initial convergence maker 
which forces the model to ’’lose” the initial state quickly and converge to 
the ’’true ocean”. The global rms-errors in the studied domain between the 
’’true ocean” and the temperature fields of the comprehensive experiment 
EX S decrease by about 50 % after 30 days and 66 % after 60 days of as- 
similation. The rms-errors of two other experiments with diagonal slices are 
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reduced by 33 %, and of all other sides of the parallelogram are reduced by 
27 % after 30 days of insertion. This shows that inserting the ” right” data 
into a "wrong” model tends to drive the ” wrong” model closer to the ’’truth” 
and to reduce global rms-errors. The good convergence also shows that the 
’’true ocean” and the model with CAT data insertion are strongly correlated. 
The rms-error curve is very close to linear in the 30 days of insertion. Then, 
it is expected that the convergence to the ’’true ocean” will be better if the 
temporal and spatial resolutions of the CAT insertion are higher. 

The temperature is spread out greater in the shallow water which may 
be the result of stronger mixing and horizontal shear and its interaction with 
topography, which increases horizontal diffusion. It is seen in all experiments 
that the temperatures of the mixed layer level at day 30 converge very well 
to the ’’true ocean” in comparison with the surface levels. This results from 
the fact that the mixing is strongest inside mixed layers. The turbulence 
activities are strongest inside mixed layers and weakest at the surface. Inside 
mixed layers both ’’top down” and ’’bottom up” processes must be active 
while only the ’’bottom up” process is active at the surface. These processes 
are important in CAT data assimilation because of insertions by vertical 
slices. 

The behavior of the global volume rms-errors is nearly insensitive to the 
change of the CAT data slice location in the more energetic regions of the 
Gulf. Although side-slices and diagonal-slices of the parallelogram have var- 
ious temperature advection, their global rms-errors distributions are almost 
the same. Their rms-errors are reduced by 27 % and 33 % after 30 days of 
assimilation. There is no the most effective slice between different local CAT 
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slices. All slices have almost the same ’’weight” (effectiveness) in improving 
the model errors on the global space scale of the Gulf of Sirte. The conver- 
gence is not proportional to the number of CAT slices in an experiment. 

The experiments also show that the internal forcing provided by the single 
slice of inserted CAT data induces an immediate correlation between the fully 
nonlinear primitive equation model and the ’’true ocean” that is not localized 
in physical space. This result is very different from that of Malanotte-Rizzoli 
and Holland who found that the single section of inserted data is too local- 
ized in physical space to induce any immediate correlation between a QG 
model and the reference ocean (MRH, 1988). The reason for our different 
conclusions is due to the data and models used and to our smaller coastal 
domain. 

Overall, a direct insertion scheme for CAT data works very well with a 
coastal multilevel, primitive equation model with complete thermodynamics, 
a second order turbulence closure, and realistic coastlines and topography. 
A comparison of this scheme with a modified Kalman filter using CAT is 
recomended for future investigation. 
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7 Figure Captions 



Table 1 Block diagram of the simulation twin experiments for studying 
coastal acoustic tomography data assimilation. 

Fig. 1 Temperature field of the ”true ocean” at the sea surface (level 1) on 
day 1. 

Fig. 2 Same as figure 1, but for level 6 (a mixed layer level). 

Fig. 3 Temperature field of the ’’true ocean” on day 1 for level 15 (the 
bottom layer). 

Fig. 4 Temperature of the ’’true ocean” at the sea surface (level 1) on day 
30. 

Fig. 5 Same as Fig. 4 but for level 6. 

Fig. 6 The ’’true ocean” at level 15 (the bottom layer) on day 30. 

Fig. 7 Temperature of the sea surface for day 1 in the experiment EXo with 
no data assimilated. 

Fig. 8 Same as Fig. 7 but for day 30. 

Fig. 9 Root-mean-square (rms) values of the differences between the exper- 
iment with the insertion of five CAT slices ( EXs ) and the experiment 
without data insertion (EXo)). 
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Fig. 10 Temperature at the sea surface on day 1 of the experiment EX t 
with a single NE-SW tomographic section. The tomographic section 
can be seen as a chain of 11 squares. 

Fig. 11 Same as Fig. 10 but for level 6 (a mixed layer level). 

Fig. 12 Same as Fig. 10 but for level 15 (bottom layer) 

Fig. 13 Temperature at the sea surface for day 30 of the experiment EXj 
with a single NE-SW tomographic section. 

Fig. 14 Same as Fig. 13 but for level 6. 

Fig. 15 Same as Fig. 13 but for level 15. 

Fig. 16 Global volume rms error versus time for a period of 30 days be- 
tween the nowcast temperature fields in EX t and the ’’true ocean” 
temperature fields. 

Fig. 17 Temperature at the sea surface (level 1) on day 1 of the experiment 
EXz with a single NW-SE tomographic section. The tomographic sec- 
tion can be seen as a chain of 11 squares. 

Fig. 18 Same as Fig. 17 but for level 6 (a mixed layer level). 

Fig. 19 Same as Fig. 17 but for level 15 (bottom layer). 

Fig. 20 Temperature at the sea surface on day 30 in EX 2 with a NW-SE 
tomographic section. 

Fig. 21 Same as Fig. 20 but for level 6. 
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Fig. 22 Same as Fig. 20 but for level 15. 



Fig. 23 Global volume rms error versus time for a period of 30 days be- 
tween the nowcast temperature fields in EX& and the ”true ocean” 
temperature fields. 

Fig. 24 Day-1 sea surface temperature of the comprehensive experiment 
EXs with five tomographic sections forming a parallelogram enclosing 
a coastal region. The parallelogram (without a southern section) can 
be seen as five chains of brightened squares. 

Fig. 25 Same as Fig. 24 but for level 6. 

Fig. 26 Same as Fig. 24 but for level 15. 

Fig. 27 Sea surface temperature on day 30 in the experiment EXs with the 
assimilation of five tomographic sections. 

Fig. 28 Same as Fig. 27 but for level 6. 

Fig. 29 Same as Fig. 27 but for level 15. 

Fig. 30 Global volume rms error versus time for a period of 30 days be- 
tween the nowcast temperature fields in EXs and the ’’true ocean” 
temperature fields. 

Fig. 31 Temperature at a mixed layer level (level 6) in the comprehensive 
experiment EXs right after 60 days of data insertion. Temperature 
of the ’’true ocean” for the same time (i.e., day 60) and same level is 
shown in the lower panel. 



Fig. 32 Global volume rms error versus time from day 31 to day 60 between 
the nowcast temperature fields in EX 3 and the ’’true ocean” tempera- 
ture fields. 



44 



Twins Experiments 
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